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ABSTRACT 


The process of selecting optimal settings for directional over-current relays 
(DOC relays) is a selection of time dial setting (TDS) and Ip (backup current), 
So that changes in the system of electrical power distribution. In this work, a 
breeder genetic algorithm (BGA) has been applied to optimal settings of DOC 
relays with multisystem D-FACTS devices. The simulation consists of two 


network operation scenarios, scenario without D-FACTS which consisting of 
coordination of DOC relays against three phase faults, and the second 


Keywords: scenarios with multi TCSC. In general, had been verified on optimal settings 
. of relays that the impacts of TCSC insertion in 33-bus distribution system on 
Breeder genetic algorithm 
DOC relays. 
D-FACTS 
Distribution system 
DOC relays 
TCSC This is an open access article under the CC BY-SA license. 
OLO 
Corresponding Author: 
Lazhar Bougouffa 


Department of Electrical Engineering 
University of El-Oued, El-Oued 39000, Algeria 
Email: bougouffa-lazhar @univ-eloued.dz 


1. INTRODUCTION 

Power system shocks are reduced by providing insulation to the faulty area as quickly as possible, 
and this is the main goal of protecting the system. The protection is operating only after an abnormal or 
insupportable condition occurs, with sufficient signal to allow their operation. Relays operate to protect the 
primary zone and also operate correctly in response two conditions outside the primary protection zone. 
Therefore, backup protection is available for the area outside its main zone in these cases. The process of 
applying and adjusting the protective relays that exceed the primary relays so that they operate as quickly as 
possible in their area, but they do not work in a faster time in their backup area. This process is known as the 
relay coordination [1]. 

In order to study the negative or positive impact on the protection system in an electric network 
when a new compensation system is introduced into the distribution system, such as the distribution 
generation resource (DGR) [2] and distributed flexible AC transmission system (D-FACTS) and improving 
the new setting of the direction over-current relays after compensation. However, with the increasing 
integration of FACTS into the distribution network, the problem of coordinating protection relays becomes 
more difficult. In recent years, researchers have provided solutions and several methods to overcome the 
effect of a series controlled FACTS device on the optimal coordination of the over-current directional relay 
problem. A new idea is introduced in order to find out how affects of the series compensation degree on the 
setting of the over-current relays. For the same problem, the authors in [3] was present a solution to the new 
setting of DOC Relays and the problem of coordination using PSO technique. The aim of this optimization is 
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to find optimal settings of TDS and Ip for relays, which reduces the operating time of the overall relays, other 
work [4] applied LP-Dual simplex algorithm to obtain the optimal TDS for a fixed Ip. L. Bougouffa et al. [5], 
the authors study the effect of DG as well as a thyristor controlled series capacitor (TCSC) on the short 
circuit level in the distribution network. 

The different optimization algorithms which were developed by researchers to investigate the 
impacts of FACTS on optimal settings and coordination of protection: Zahra M, Mostafa J and Mehdi Gh. [6] 
Modified adaptive particle swarm optimization (MAPSO), N. Mancer, B. Mahdad, K. Srairi et al. [7] Particle 
swarm optimization (PSO) algorithm based on time-varying acceleration coefficients (PSO-TVAC), E. 
Dehghanpour, H. K. Karegar, R. Kheirollahi et al. [8] Cuckoo-linear optimization algorithm, R. Benabid, M. 
Zellagui, A. Chaghi et al. [9] Differential evolution algorithm, Behzad A. M. A., Masoud E. T., Reza R. et al. 
[10] An adaptive cumulative sum based method. 

In this work, breeder genetic algorithm (BGA) is proposed to select the optimal values of TDS for a 
fixed value of Ip and present a solution for the coordination problem of DOC relays in presence of D-FACTS 
devises. In this paper, the algorithm is applied to IEEE 33-bus system with multi-TCSC which is modeled 
and simulated to prove the efficiency of the proposed BGA. Moreover, the obtained results when using this 
algorithm are compared without and with TCSC. 


2. COORDINATION OF RELAYS PROBLEM 

Typically, the inverse time over-current relay consists of two elements, an instantaneous unit, and a 
time dependent unit. A time delay unit is used for the current, which is less than the instantaneous current 
setting, but exceeds the normal flow due to a fault. This unit operates when malfunction with a time delay. 
The two decisive factors are TDS and IP. The IP is the minimum current value for which the relay operates 
and the TDS defines the tripping time [11], [12]. In (1), the standard mathematical formula for modeling the 
operating time of relays according to IEC 60255 [13]: 


t = TDS x 


Where, a and B are constants. The values of these constants are given in [14], [15], Ip is the current in 
abnormal operation at the primary of a current transformer, CTr is the current transformer ratio, DOC relays 
will operate rapid after sensing a high current in abnormal operation because the tripping time of relay is 
inversely proportional to the fault current [16]. 


2.1. Objective function 

The TDS setting and the IP setting are two parameters that define the optimal coordination problem. 
Conventionally, in coordination studies, the objective function is formed as the sum of the operating times of 
all the primary relays [17]. In (2) represent the minimization of total operation time [6]: 


= 0.14 
i=1 eee 1 
Cl. T) 7 


Where OF is the objective function in sec, each relay would have 0.2 seconds as minimum operation time 
[18]. 


2.2. Bonds on operation time of relays 

The aim of this study is to ensure that the over-current relays reduce faults in a very fast time, while 
maintaining coordination between the backup protection zone and the primary protection zone relays. 
Coordination is maintained when relays in primary protection zone operate faster than relays in backup 
protection zone with a difference time equal to or greater than the specified coordination interval. Three 
constraints are considered for the minimization problem [12]; the first constraint is CTI; 


T nii -up ` T PR <CTl (3) 
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The CTI should be the minimum margin between the pair of relays primary and backup [19, 20]. The margin 
of safety is set between 0.2s and 0.5s depending on the circuit breaker speed and the type of relay [21]. 


2.3. Bonds on TDS and Ip for each relay 
The main decision variable in the traditional approach to coordination of protection is the TDS. In 
(4) the TDS value for the relay i is limited between the TDSi”” and the TDSi”®™ [22], [23]. 


TDS?" STDS ; STDS ;” (4) 


|e BS hag (5) 


Where: TDSiE [0.05, 1.5] and pE [0.25:0.25: 2.5]. 


3. GOAL OF STUDY AND IMPACTS OF D-FACTS 

Directional over-current protection is used to protect the network, in abnormal operations in the 
distribution network. The presence of D-FACTS disturbs the performance of distribution network protection. 
However, the total fault current seen by the relays will be changed when there are D-FACTS. Based on the 
TCSC operation in inductive and capacitive modes, the Fault current might increase or decrease [24], [25]. 
The intensity of the current measured by the relays changed where TCSC is connected to the distribution 
network. Also, DOC relays are incapable to detect the fault current because TCSC operate as a fault courant 
limiter. 


4. BREEDER GENETIC ALGORITHM 
The BGA is a new evolutionary algorithm and similar to genetic algorithm (GA). However, BGA is 
based on artificial selection unlike GA which is based on natural selection. 


4.1. Selection 

The BGA uses truncation selection method. To perform this method, trial solutions are evaluated 
and ranked according to their fitness. The top T% of the best individuals of the trial solutions is selected and 
goes through recombination and mutation to create the next generation. The rest of the individuals are 
discarded. The next generation of the best test solution, which is directly inserted in the next generation and 
the rest of the solutions is created by recombination and mutation. Recombination is similar to crossover in 
GAs. The process is repeated until an optimal solution is obtained or the maximum number of iterations has 
been reached. BGA is very versatile and easy to program [26]. 


4.2. Recombination 
In volume recombination, a random vector r equal in length to the parents’ length is generated. The 
child zi is produce by the expression [26], [27]. 


7a A); (6) 


Where xi and yi are two parents. 
In line recombination, a single uniformly random number r is generated between O and 1, and the child is 
obtained by the expression [27]. 


Ff a oe an oe ea 8 8 (7) 


4.3. Mutation 

One problem that has been of concern in GAs is premature convergence, where the search might 
converge on local optima rather than the desired global one. This has been minimized by preserving the 
diversity of the population by adding small injection of randomness or mutation [27]. This is achieved by 
adding a small vector of normally-distributed zero mean random numbers to each child before inserting it 
into the population. The standard deviation r of the vector is very critical, as small r might lead to premature 
convergence or big r might impair the search and reduce its ability to converge optimally. Therefore, it’s 
better to use an adaptive approach were-by the rate is modified during the course of the search. The 
population is divided into two halves X and Y. A mutation rate of 2r is applied to X while a mutation of Yr to 
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Y. The mutation r is adjusted depending on the population (X or Y) that is producing better and fitter 
solutions on average. If X individuals are fitter that the mutation rate, r is increased by 10%, while if Y is 
fitter that the mutation rate, r is reduced by similar amount. 


5. DETERMINE THE RELAYS TDS 

The TDS adjusts the temporization before the relays operates whenever the current measured 
reaches the value equal to or greater than current setting of the relay [14]. In order to obtain appropriate 
protection and coordination, the proposed algorithm is illustrated in the Figure 1. 


Initialization the data of network and TCSC 





Calculate the fault with and without TCSC 
od Compute TH primary, TH Backup and CTI 


e Evaluate the fitness of each TDS 


population using the objective function 


No 
Termination 


Yes 


e Generate initial population of TDS 


e Round each TDS to the available 
integer 





e Compute T" primary, T“ Backup and CTI 
e Evaluate the fitness of each TDS 
population using the objective function Display Values of TDS corresponding to 


lowest objective function value 


e Check Constraints 





i=i+] | Select TDS Values through truncation 
selection for recombination 





Create new TDS values 


Figure 1. Flow chart for application of BGA for relay coordination with TCSC 


5.1. Power system model used for studies 

The network diagram used for studies is given in Figure 2 [3], [28]. BGA was applied on the 
12.66 kV, 100 MVA distribution systems which consist of 33 buses, 32 directional over-current relays and 
four TCSC. Before the BGA is applied to the problem of coordination with TCSC, the following has to be 
done: 
— Determine the pairs of relays for coordination 
— Determine the fault current for each relay pair with and without TCSC. 
— Determine the TDS for each relay for a fixed the pickup current to 1. 
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Figure 2. System diagram of power system model used for study 


The validity of BGA is improved and tested on two scenarios. The first study consists of calculate 
the optimal settings of DOC relays without TCSC. Second case study consists of determine the optimal 
settings of DOC relays with the presence of multi-TCSC. The parameters of the TCSC are L=9.60 mH, 
C=212.2 uF. The rated value of TCSC is a function of the reactance where the TCSC is installed and 
expressed as: 


Xiine = Xij + Xresc (8) 
Where: 
Xrosc = Krcesc: Aline (9) 


Where, Xjineis the overall line reactance with TCSC installation. x;;Is the reactance between bus-i and jJ. 
Xrcsc 18 the reactance of TCSC and Kyrcsc is the coefficient which represents the compensation level of 
TCSC. The working range of reactance of TCSC is fixed between -0.7 (capacitive) Xine and 0.2 (inductive) 
Xjine [29]. The TCSC is installed in the middle of the line; the three-phase fault equals the direct component 
is given by: 


V + Vrese 
l = 








Z. Z. (10) 
Zij-1 + ( 5 +) + Xrcsc-1 + ( 5 +) 


The opposite and homopolar components are equals to zero: 
The current of three phases are given by: 
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(V; + Vrese) 


pae a 
AA Za (11) 
Zij-1 + ( 5 +) + Xrcsc-1 + ( 5 +) 








eos (V; + Vrese) (12) 


ae Te, 
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Zij-1 + (==) + Xrcsc-1 + ( 5 











Where, Z;;_,1s the direct component between bus-i and j. Vrcsc is the voltage of TCSC 


6. RESULTS AND INTERPRETATIONS 

The results are provided to demonstrate the performance of the relays to isolate faults as quickly as 
possible for primary and back-up protection. The performance is evaluate by looking at the coordination time 
interval, operating time for primary zone fault and operating time for backup zone fault. 

We have selected two angles for each mode (inductive and capacitive). Finally, the performance of 
settings for relay pairs is set in Table 1. The fitness value for the BGA algorithm drops rapidly in the 
beginning of the search from the initial value of 700 seconds to 78 seconds where it begins to settle to a 
fitness value in 10 iterations. It settles between 78 and 70 seconds for about 20 iterations. The algorithm 
maintains diversity and converges steadily until it reaches a fitness value that is within 7%. 


Table 1. Optimal setting of TDS using the BGA 
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With TCSC 
mig Without inductive capacitive 
TCSC 90° 135° 136° 180° 

TDS! 0.5200 0.5000 0.9700 0.1800 0.6300 
TDS’ 1.2400 1.3400 1.6000 0.3400 0.5000 
TDS? 1.0600 0.0900 1.4400 1.1600 0.3800 
TDS* 0.9200 0.1400 1.1400 1.1800 0.5500 
TDS” 1.1300 0.8200 1.4600 0.4700 0.5200 
TDS 1.1500 0.1000 1.0100 0.8200 0.3700 
TDS” 1.1200 0.5200 1.3800 1.0000 0.4700 
TDS’ 1.0800 0.1400 1.2000 0.5800 0.1800 
TDS? 1.0800 1.1100 1.5000 0.8100 0.1700 
TDS”? 1.0600 1.0900 1.0000 1.0100 0.5300 
TDS” 1.0300 0.1800 1.3100 1.0800 0.5400 
TDS” 1.2000 0.0900 1.1900 0.8600 0.1100 
TDS"? 1.1100 0.1000 1.0700 0.5400 0.1900 
TDS" 1.0400 0.1400 0.9800 1.2000 0.2300 
TDS” 1.1300 0.3000 1.0800 0.4600 0.1200 
TDS’? 1.1800 0.1900 1.0700 0.4100 0.3700 
TDS?” 1.0500 0.1800 1.2800 0.8500 0.1100 
TDS!’ 1.1600 0.1700 1.4100 0.7300 0.3500 
TDS”? 0.9700 0.6400 1.2500 1.1200 0.1200 
TDS” 1.2400 0.1900 1.2900 0.8400 0.1300 
TDS” 1.1800 0.2100 1.0400 0.3100 0.4700 
TDS” 1.1000 0.1000 1.1600 1.0100 0.2000 
TDS” 1.1400 0.3700 1.5700 0.4500 0.4000 
TDS” 1.1300 0.6800 1.3700 1.1200 0.3000 
TDS” 0.7200 0.1900 1.1500 0.5200 0.4700 
TDS *° 1.0300 0.1900 1.1300 1.1400 0.3000 
TDS” 1.0900 0.1900 1.1000 0.7800 0.1100 
TDS” 1.1100 0.1700 1.3500 1.1400 0.4400 
TDS” 1.1300 1.0800 1.0200 1.0000 0.5000 
TDS? 1.1600 0.4200 1.1800 1.1800 0.1800 
TDS?! 0.8900 0.3800 1.1700 0.1700 0.3500 
TDS”? 1.2400 0.3300 1.1900 1.0700 0.5000 
Obj.F 75.9454 75.8682 76.8308 78.1341 77.3951 
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We see only the variations of the short-circuit current as a function of the reactance injected by the 
TCSC (ie substantially the operating mode of the TCSC in the system) is minimal in the inductive mode and 
maximum in the capacitive mode. According to the results obtained by the BGA method, it can be concluded 
that the TDS setting of all relays is less than 1.2 and greater than 0.05, hence a satisfied standard. The main 
relay for the selected pair is Relay 8 and the backup relay is Relay 7. For this relay pair, three phase fault was 
simulated in bus-barre 9. This fault is in the zone of protection for relay 8 and in the backup zone of 
protection for Relay 7. For the proposed BGA, the TDS of relay 8 is 1.0800 seconds and Relay 7 is 1.1200 
seconds without TCSC located at the branch between bus-barre 8 and 9. In the presence of the TCSC and its 
modes of operations, the results of the optimal settings of the pair of relays is clear in the Table 1 to guarding 
the coordination time interval of 0.3 seconds. 


7. CONCLUSION 

This study allowed us to design a procedure for optimizing the parameters of DOC relays in the 
presence of multiple TCSCs, based on BGA. This is performed to obtain the optimal setting values of TDS 
for each relay in IEEE 33-bus distribution system. The optimal TDS values achieved the required minimum 
operating time, which is 0.2 s for different TDS values. BGA show high capability to solve the coordination 
problem and provide better results compared to conventional methods. 

We were able to achieve the objectives set out by obtaining efficient and robust relays. There are 
still many avenues to explore, but it seems to us a priority to continue a more in-depth study on the 
application of these techniques for simultaneous multi-objective optimization of D-FACTS and their 
locations by considering criteria such as minimization of response time, and the cost of installing D-FACTS. 
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